Quantum plasma effects in the classical regime 
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Abstract 

For quantum effects to be significant in plasmas it is often assumed tliat tlie temperature over 
density ratio must be small. In this paper we challenge this assumption by considering the con- 
tribution to the dynamics from the electron spin properties. As a starting point we consider a 
multicomponent plasma model, where electrons with spin up and spin down are regarded as dif- 
ferent fluids. By studying the propagation of Alfven wave solitons we demonstrate that quantum 
effects can survive in a relatively high-temperature plasma. The consequences of our results are 
discussed. 

PACS numbers: 52.27.-h, 52.27.Gr, 67.57.Lm 
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Recently, several studies of quantum plasmas have appeared in the literature 
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Quantum plasma effects can also be seen in scattering experiments with solid 



density targets [24|. The usual regime where quantum effects are important involves dense 
low-temperature plasmas, where either the Fermi pressure is comparable to the thermal 
pressure or the thermal de Broglie wavelength times the plasma frequenc y is cornparable 
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to the thermal velocity. In recent studies of spin effects in plasmas |8|,l9|,ll0|,lll|,ll2|, the 
condition for quantum effects to be important has been found to be somewhat different from 
the case of non-spin quantum plasmas , but also here a high temperature 

tend to make quantum effects small. 

In the present Letter we study a weakly collisional high temperature plasma. In partic- 
ular, we focus on the case where the temperature over magnetic field ratio is sufficiently 
high to make spins randomly oriented at thermodynamic equilibrium. Within the one-fluid 
model, such a condition tends to make the macroscopic spin effects negligible 
However, here we study a two-electron fluid model, where the different electron popula- 
tions are defined by their spin relative to the magnetic field. Evaluating this model for the 
particular case of Alfven waves propagating along the external magnetic field, it is found 
that linearly the predictions agree with the one-fluid spin model. Nonlinearly, however, the 
induced density fluctuations of the spin-ponderomotive force is significantly different for the 
two-spin populations. As a consequence, the self-nonlinearity of the Alfven waves gets a large 
contribution from the spin effects, even for a high-temperature plasma. In general, the con- 
clusion is that spin effects cannot be neglected even in moderate-density high-temperature 
plasmas that normally are regarded as perfectly classical. 



Neglecting spin-spin interactions, the equations of motion are 



dtUs + V ■ {us^s) = 0, 



0, y, 0, 
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where Us and v<j are the density and velocity of species s, s = i, +, — enumerates the plasma 
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particle species and ± denotes the two types of electrons, 



msUs [dt + Vs ■ V) = QsHs (E + X B) - c^Vus 




(2) 



and 



{dt + V, ■ V) S, 
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where is the spin of species s, qs is the charge of species s, ps = PsiTg] Us) is the pressure 
of species s and Cg = {dps/ dnsY^"^ is the sound speed of species s (where we have assumed 
an isothermal plasma) containing also contributions due to the Fermi pressure, /i^ is the 
magnetic moment of species s, and — yU± = fiB = eh/2me is the Bohr magneton, e is the 
magnitude of the electron charge, h is Planck's constant, me is the electron rest mass, and 
c is the speed of light. We note that Einstein's summation convention has been used in (2). 
In what follows we will, due to the large ion inertia, neglect the quantum corrections to the 
ion motion. 

For temporal variations of the magnetic field faster than the inverse electron cyclotron 
frequency, spin flips can be induced. Furthermore, particle collisions can also reverse the 
spin. Thus, to make sure that spin reversal does not occur, we consider dynamics on a 
time scale shorter than the inverse collision frequency, but longer than the inverse cyclotron 
frequency. For this case, we can replace the spin evolution equation ([3]) with the relation 
S-t = =p(/;,/2)B, for electrons with spin up and down relative to the external magnetic field. 

The coupling between the quantum plasma species is mediated by the electromagnetic 
field. The magnetizations due to the different spin sources are M± = —2fiBiT'±^±/fi = 
±HBn±B. Ampere's law then takes the form 





where the free current is denoted j. Moreover 




(5) 



The system is closed by Faraday's law 



V X E 



dfB. 



(6) 
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In previous works [8|, [o], Q, [ill, [1^, electrons have been treated as a single population, 
with a single macroscopic velocity v and spin vector S. As argued above, for fast dynamics 
in an approximately collisionless plasma, this is not appropriate, as the populations with 
positive and negative spins along the magnetic field will not change spins on the short time 
scales considered, and as seen in Eq. ([2]) the two populations are described by separate 
evolution equations. If we describe the plasmas by a single electron population, with a 
background spin distribution close to thermodynamic equilibrium, the spin effects are limited 
to a certain extent whenever /i^Bo/A^fiTe <^ 1. This is due to the thermodynamic Brillouin 
distribution for spins oc taiih.{fiBBo/kBTe) describing the macroscopic net effect of the spin 
orientation. Thus within the single electron fluid model, we need low temperatures or 
very strong magnetic fields for spin effects to be important. By contrast, within the two- 
fluid electron model, spin effects may be of importance also in a weakly magnetized high- 
temperature plasma, as will be demonstrated below. 

As an example, we consider the nonlinear response to a low-frequency electromagnetic 
Alfven wave pulse propagating parallel to an external magnetic field. In linear ideal mag- 
netohydrodynamic (MHD) theory, the magnetic field perturbation thus propagates along 
the external magnetic field Bq = Bqz with the Alfven velocity ca = (-Bg/z^oPo)^^^) where 
Po is the unperturbed mass density. Since linearly the Alfven wave has no density per- 
turbation, the quantum effects of the Bohm-de Broglie potential and the Fermi pressure 
do not change this result. In what follows, such quantum effects will be omitted. As 
found in within a single fluid spin model, the Alfven velocity is decreased by a factor 
1 + {hijp^/2mic'^ujie')ta.nh.{jj,BB()/kBTe) due to the spin where Upe is the electron plasma 
frequency, rrii the ion mass, cuce'' = e/io-ffo/'"^e is the electron cyclotron frequency due to the 
magnetic field fioHo = Bq — /ioM], which is the field with external sources only, i.e. the 
contribution from the spins are excluded (here Mq is the unperturbed magnetization due to 
the spin sources). For pbBq <^ kBTe the correction factor for the Alfven velocity is close 
to unity, and the approximation fioHo ^ i?o is a good one. This is the case that will be 
considered below, and the spin corrections to the linear Alfven velocity will therefore be 



omitted in what follows. Furtermore, the envelope of a weakly modu 



ated Alfven wave will 



propagate with a group velocity Vg ~ ca for frequencies u <ti uj^ |25|, where Ud is the ion 
cyclotron frequency. 

The ponderomotive force of this envelope will drive low-frequency longitudinal pertur- 
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bations (denoted by superscript 'If in what follows) that are second order in an amplitude 
expansion, and to leading order depend on a single coordinate ( = z — CAt- Thus the dynam- 
ics is considered to be slow in a system comoving with the group velocity. Neglecting spin, 
for a slowly varying magnetic field perturbation of the form B = 5(2;, t) exp[i(A;z — ujt)]e 
+C.C. (where e is a unit vector perpendicular to z and c.c. denotes complex conjugate), the 
low-frequency MHD momentum balance can be written 

p^dtv'l = -d, [\B\^ //io + (M + kBT,)n''] , (7) 

where index i denotes ions. For simplicity we assume a weak magnetic field, Bq ^ 
[jJoPo{kBTi + kBTf.) /'mi\^/'^ . Relating the low-frequency perturbations of the density and 
velocity using Eq. ([T]), the left hand side of Eq. ([7]) is then found to be negligible, and the 
density depletion is given by 



\B 
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Moreover, since we have charge neutrality within the MHD approximation, we have here 
neglected the index i on the density, since the total electron density will be the same for a 
proton-electron plasma. Next, we add the spin terms in our model. Again neglecting the 
ion inertia, the MHD low-frequency force balance Eq. ([7]) is replaced by 

Fp+ + Fp_ = d-, [\B\^Ip^ + {kBTi + kBT,)n''] (9) 

where the ponderomotive force contributions Fp-t are low-frequency perturbations due to 
the terms 2fisnsS^'VBa/h in the electron momentum equations. Including the spin vectors 
component in the direction of the perturbed magnetic field, the sum of these spin force 
contributions can be written 

no uuBn eh d , , 

Fp+ + Fp_ ^ LX 10 

where we have used that the unperturbed density difference (no+— no-) of the two spin popu- 
lations in thermodynamical equilibrium is proportional to ta.nh.{fiBBo/kBTf.) ~ pBBo/kBTf.. 
Thus the net effect of the spin ponderomotive force on the ion density as well as the total 
electron and ion density is very small in the regime fiBBo/kBTe <^ 1, similarly as we would 
have for a single electron spin model, and consequently Eq. ([8]) is a valid approximation. 

The interesting difference between different fluid models comes when we analyze the den- 
sity perturbations of the two electron populations separately. The low-frequency momentum 



balance equation for the different electron species are 
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where we have introduced the electrostatic low-frequency potential (this potential does 
not appear in the overall momentum balance, as the plasma is quasi-neutral). By adding 
the ± parts of Eq. (fTTj) and integrating we obtain 



to first order in the expansion parameter fiBBo/kBTe. Due to the small factor fiBBo/kBTg 
in front of the last term, the spin effects on the total electron population are small, in 
agreement with Eq. (fTOj) . again justifying the omission of spin effects in Eq. (|8]). However, 
solving instead for the density difference between the two electron populations we find 

(13) 
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The importance of the density difference displayed in (fT3!) appears when the nonlinear 
self-interaction of the Alfven waves is studied. The momentum equation contains the term 
(j X B/n)'^\ where nl denotes the nonlinear part and j is determined from Eq. (jlj) with 
the displacement current neglected. Omitting the terms of higher order in the expansion 
parameter fiBBo/kBTe we find 



j X B 
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where the last term represent the two-fluid electron spin contribution. This results in a 
corresponding spin- modification of the self-nonlinearity of the Alfven waves. Including 
weakly dispersive effects due to the Hall current, parallel propagating Alfven waves are 
described by the so called derivative nonlinear Schrodinger (DNLS) equation 26!]. For a 
quasi-monochromatic wave as considered here, the DNLS equation reduces to the usual 
nonlinear Schrodinger (NLS) equation [23] of the form 

(15) 



2 ' Bi 







Here the group dispersion v' = dvg/dk is the group dispersion, ( 



Vgt is the comoving 



coordinate and Vg is the group velocity. These quantities are determined from the Alfven 



wave dispersion relation, which reads a;^ = k'^c\{l±kcA/uJci), when weakly dispersive effects 



due to the Hall current is included 



25|. The upper (lower) sign corresponds to right (left) 



hand circular polarization. The nonlinear coefficient is Q = Qd^ — i'^f^BBo/mic\y], where 
the classical coefficient is Qc = kc\/4:{c\ — cj.) ~ —kc\/4:cj.. The NLS equation has been 
studied extensively |[28|], and as is well-known it admit soliton solutions in ID, and can 
describe nonlinear self-focusing followed by collapse in higher dimensions. Furthermore, the 
evolution depends crucially on the sign of the nonlinear coefficent, which may change due 
to the spin effects. However, our main concern in this context is not the evolution of the 
Alfven waves, which were chosen just as an illustration. The fact that interests us here is 
that spin can modify the dynamics even when the spins are almost randomly distributed 
due to a moderately high temperature (i.e. when /xbSo/^b^g ^ !)• The approximately 
random distribution of spins is shown in the dispersion relation of the linear wave modes, 
which are more or less unaffected by the spins since linearly the total spin contribution 
on the electrons cancel to leading order. Nonlinearly, however, the consequences of the 
different density fluctuations induced in the spin-up and spin-down populations are seen. 
The unique feature of this quantum effect is that is survives even for a high temperature. 
By contrast, well-known quantum plasma effects like the Fermi pressure, and the Bohm-de 
Broglie potential becomes insignificant for high temperatures. This is also true for single 
fluid spin effects js, 9, 
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An illustration of the regimes where the different quantum plasma effects become signifi- 
cant is provided in Fig. 1. In particular we note that the two-fluid nonlinear spin effects are 
important for high plasma densities and/or a weak (external) magnetic fields. For compar- 
ison, both the Fermi pressure and the Bohm-de Broglie potential need a low-temperature 
or a very high density to be significant. Single-fluid spin effects can also be significant in 
this regime, or in the regime of ultra-strong magnetic fields that can occur in astrophysical 
applications. Especially interesting is that two-fluid nonlinear spin effects can be significant 
in a high-temperature regime that is normally perceived as classic. While this obviously is 
an intriguing result, a word of caution is needed. Although our results clearly show that 
spin effects can be important when fiBBo/niiC^ approaches unity, we note that in a number 
of applications, spin effects of the kind discussed here can be suppressed even if ^BBo/mic\ 
is large. These include: 

— Systems where the dynamics is dominated by compressional effects. In such problems 
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the thermal pressure force dominates over the spin force of the electrons, and spin effects 
are suppressed unless fiBBo/kBTe 1; 

— High density collisional plasmas. In such cases the spin orientation changes rapidly 
and different spin populations cannot be sustained; 

— Systems where the spin forces are negligible as compared to the electrostatic force. 

Furthermore, the strong magnetic fields needed for confinement tend to make the two- 
fluid spin effects studied here negligible for magnetically confined plasmas. Nevertheless, 
even with the above cases excluded, our discussion shows that there is a large range of plasma 
problems that traditionally have been dealt with using purely classical plasma equations, but 
where the electron spin properties give a significant contribution to the dynamics. In general 
the mechanism can be summarized as follows: For a weakly magnetized initially homoge- 
neous plasma, the spin-up and down populations are (approximately) equal. However, when 
an electromagnetic perturbation enters the system, the spin ponderomotive force separates 
the two populations, which in turn modifies the magnetic field since spin-magnetization no 
longer cancels. From then on, a two-fluid electron model is needed. In particular, in the 
region of aligned electron spins, the original magnetic field will be enhanced, and hence 
mechanisms of this type can play the role of a magnetic dynamo. Suitable plasma condi- 
tions for nonlinear two-fluid spin effects to be important may be found in weakly magnetized 
inertially confined plasmas, near atmospheric pressure plasma dischages, weakly magnetized 
regions of the sun's convection zone, as well as in astrophysical and cosmological plasmas. 
Exploration of these vast range of problems remain a major research project. 
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